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ABSTRACT 



The determination of overpressure yield from internal 
blast has been the subject of prior investigations. Models 
have been developed to predict explosion thermodynamic param- 
eters, based on equations of state, mass and energy conserva- 
tion, and rate laws, utilizing high speed digital computers. 
Typically, thermodynamic equilibrium was assumed, an equilib- 
rivim temperature was determined, ideal gas law behavior was 
ascribed to product gases, and a pressure was computed. 

Reports of the results of these studies have typically con- 
tained a disclaimer that radiation effects have been ignored. 

In this investigation, a blackbody radiator model was developed 
and, utilizing the parameters reported in the literature for a 
number of explosives of common interest, a determination was 
made of the effect of thermal radiation as an attenuation 
mechanism on explosive energy and the existence of thermody- 
namic equilibrium. Energy losses of less than .5% were calcu- 
lated and thermodynamic equilibrium was demonstrated through 
statistical considerations. 
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I . INTRODUCTION 



Internal blast is the term applied to the concussion 
effect resulting from highly exothermic reactions occurring 
in a confined space. Examples of internal blast are dust 
explosions in grain elevators and coal mines; gas pocket 
explosions in buildings and underground systems such as 
sewers and mines; and the detonation of high explosives, 
either accidently or through proper fuse function, in a 
charge hole in a mine ore face, or in a shipboard ordnance 
stores magazine or other compartment. Substantial effort 
has been devoted to the study of explosive phenomena during 
the past century; internal blast investigation has been 
particularly intense during the past two decades. 

Various computer codes have been developed to predict 
the thermodynamic properties of the detonation products 
generated by explosive materials. Typically, these codes 
have been used to determine a product composition and condi- 
tion based on chemical reaction rate equations and other 
thermodynamic and conservation considerations applied to a 
system of equations of state, such as the STRETCH code, which 
is based on the Becker-Kistiakowsky-Wilson (B-K-W) equations 
of state, reported by Mader of the Los Alamos Scientific 
Laboratory [Ref. 1) . 

Other work has been very specifically directed to the 
determination of internal blast overpressure yield. Previous 
investigations were conducted by G. F. Kinney, R. A. Reinhardt 



6 



and K. J. Graham of the Naval Postgraduate School, working 
in conjunction with personnel from the Naval Weapons Center, 
China Lake, California, on determination of peak overpressure 
yield of various explosives and fuels, with and without the 
incorporation of blast-augmenting reactive metals [Ref. 2, 

3 and 4] . The technique for determining internal blast 
overpressure used by Kinney, Sewell and Graham and by 
Reinhardt is to solve a system of simultaneous equilibrixam 
equations to determine a product composition, based on a 
selected trial temperature. The sum of the internal energies 
of each product gas species at this trial temperature is 
then computed, using a polynomial form for c^ for each 
species. The value of the sum is compared to the previously 
determined value of the maximum energy available for release 
in the explosion. The value of the trial temperature is 
varied until energy balance is achieved. This temperature is 
an equilibrium flame temperature. Ideal gas law behavior is 
then ascribed to the mass of product gases, and a peak over- 
pressure is computed. 

A second major effort was undertaken at the Naval Ordnance 
Laboratory, White Oak, Maryland, with a focus on the produc- 
tion of a computer program capable of predicting damage 
mechanism due to internal blast [Ref. 5] . A common disclaimer 
made in the reports of these investigations has been that the 
effects of thermal radiation as an energy loss mechanism can 
be ignored. Additionally, an assumption of thermodynamic 
equilibrium was made in the determination of flame temperature. 
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The '.validation of the assumption of thermodynamic equilibrium 
establishment in the detonation reaction zone and the magni- 
tude of the thermal radiation losses are the subject of this 
investigation. Examination of open literature reports of 
investigations conducted on lum.inosity of explosions was made 
to determine if thermal radiation energetics in conventional 
high explosive detonations had been quantified. 

Inconsistencies of values which various sources cite for 
the same parameters is noted in examining the literature 
reporting explosives experimental results. A remark by 
Rudlin [Ref. 6] to the effect that statistical deviations 
within explosives data sets ostensibly reporting replicate 
measurements average about 10 to 15%, and that agreement to 
within 5% was unusually good, provides a basis for the order 
of magnitude and two digit precision used in the model 
development for this investigation. The best experimental 
or computational values available are used, the best being 
defined as having been cited in more than one source, or 
being internally consistent with other data cited in more 
than one source. 

A model of a detonating explosive is developed to explore 
the existence of thermodynamic equilibrium and to quantify 
expected radiant exitance from the detonation. This model 
is established as a blackbody, with size and temperature 
characteristics determined by the detonation parameters of 
the explosives considered, as found in the literature. Util- 
izing these parameters, a determination is made of the effects 
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of thermal radiation as an attenuation mechanism on the explo- 
sive energy available for generation of overpressure, and of 
the existence of thermodynamic equilibrium. 
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II. DETONATION THEORY 



A brief discussion of the theory of detonation and terms 
used in detonation literature is considered appropriate prior 
to discussion of the specific model developed. Condensed 
explosives are solids which consist of materials or mixtures 
of materials which are capable of undergoing highly exothermic 
chemical reactions at extremely fast rates / producing gaseous 
and, usually, solid reaction products at high pressures and 
temperatures. "CHNO" explosives are organic materials whose 
molecules contain atoms of carbon, hydrogen, nitrogen and 
oxygen. These complex molecules undergo a decomposition reac- 
tion, followed immediately by oxidation-reduction reactions 
which produce low molecular weight detonation products, such 
as H2O, CO, CO2, N2 and C. The exact quantity of each reaction 
product, per mole of explosive, is dependent on such factors 
as the initial composition and density of the explosive, the 
fuel-to-air ratio and initial conditions of the environment. 
Detonation of an explosive results from rates of reaction 
which are so rapid that the heat evolved in the exothermic 
decomposition and subsequent oxidation-reduction processes 
cannot be carried away before there is an energy build-up of 
such magnitude that a shock front results. This shock front 
propagates through the explosive material; it is called the 
detonation wave, and propagates at a characteristic velocity 
of detonation, D. 
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The detonation shock front separates the unreacted explo- 
sive material from the material undergoing reaction. One 
description (the model adopted for this investigation) considers 
three areas in an explosive undergoing detonation, separated by 
two planes, as depicted in Figure 1. In this description, the 
shock wave initiates the decomposition reaction of the complex 
"CHNO" molecule. Decomposition reactions are completed and 
oxidation-reduction reactions proceed to completion in the 
reaction zone. When the "C-J plane" is reached, reaction 
products are at thermodynamic equilibrium, and are of composi- 
tion and condition predictable by appropriate equations of 
state. 

Part of the net energy in the reaction zone pumps the 
detonation wave to sustain it as it propagates through the 
explosive. The rest of the net reaction energy is taken up 
as internal energy of the reaction products . When the detona- 
tion wave reaches the surface of the explosive, it shocks the 
surrounding medium, setting up the initial shock wave in the 
medium, and establishes a transient reflected shock wave 
back into the material, which results in a compression. Thus, 
the net energy released in the detonation is partitioned, with 
some going into the surrounding medium as the initial shock- 
wave, and the rest being retained in the explosive material. 

The values of specific heat of detonation, q, reported 
in the literature are often based on constant volume bomb 
calorimetric measurements. Values of heat of detonation 
experimentally determined in calorimetric devices are lower 
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than the maximum heat of detonation calculated for the enthalpy 
of the reactions transforming the explosive to its detonation 
products [Ref. 8 ]. The bomb calorimeter values are lower 
than calculated values because the products formed in the 
bomb differ from those in the Chapman- Jouguet (C-J) detonation 
state. The final products achieved in the bomb calorimeter 
are not the C-J products but, rather, reflect approximate 
equilibrium at low pressures. The relative amounts of H2O, 

CO2 and solid carbon formed in the reactions, as compared to 
CO and H2f are determined by pressure dependent equilibrixam 
equations. The formation of H2O, CO2 and solid carbon are 
preferred over the formation of CO and H2 at the very high 
pressure of a C-J detonation state. The detonation products 
analyzed from bomb calorimeter tests, even when confined, 
show much higher concentrations of CO and H2 than predicted, 
at the expense of H2O, CO2 and solid carbon [Ref. 8 ]. The 
enthalpy of the C-J detonation products is more negative than 
that for the bomb products. Since the object of this investi- 
gation is to determine the thermal radiant loss from the more 
energetic detonation state, calculated values are used [Ref. 

7 ] , where applicable. 

A brief review of the pertinent points of the theories 
which have been advanced to explain explosive behavior pro- 
vides definition to the terminology and concepts presented 
above. Reference 8 , chapter 6 , contains a satisfactory devel- 
opment of the elementary theory of the plane wave detonation. 
Extension of the concepts presented there to a spherical wave 
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is straightforward [Ref. 8, ch. 9]. The detonation wave, in 
the framework of elementary thermodynamic-hydrodynamic theory, 
is asstimed to be a shock wave, treated as a mathematical dis- 
continuity, with full application of the general conservation 
laws which apply across any such discontinuity. As such, then, 
the detonation wave obeys the Rankine-Hugoniot equation for a 
shock wave; 

®1 ■ ®o = ‘■‘Pi * Po> ‘■'o ■ ■'l' 

where the subscript o refers to the unreacted material and 
1 to the reaction products; e is the specific internal energy 
function, p is the pressure and v is the specific volume. 

This equation will be referred to as the Hugoniot. If the 
products are assumed to be in thermodynamic equilibrium, then 
their thermodynamic state is specified as 

®1 ^ 

where e(p,v) is some known function of the state variables. 

In terms of specific enthalpy, 
h = e + pv, 

the Hugoniot can be written in an alternative form 

hi - ho = (Pi - Po> (Vi + Vo> . 

A plot of the Hugoniot in the plane provides a frame- 

work for definition of the Chapman-Jougtet (c-J) point (see 
Figure 2) . 

Consider the point on the Hugoniot curve for which v^^ 
is equal to v^, corresponding to conversion of the explosive 
to its products at constant volume. The value of pressure 
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PIGURE 2 
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is the minimum 



for this volume on the Hugoniot, p^, 
value of pressure for detonation, and defines the lower 
point detonation limb of the Hugoniot. Mass 

and momentum conservation conditions yield the relationship 



2 o2 

- Po ° 



(Pl - P]^) 
<''1 - ''o* 



which describes a straight line in the P^ ” plane; it is 
called a Rayleigh line, and is the chord connecting some 
point (p^,v^) to the point (p^/V^) . The propagation velocity, 
D, is proportional to the square root of the negative slope 
of the Rayleigh line. One Rayleigh line passes through 
(p^,v^) and is tangent to the Hugoniot on the detonation limb. 
The point of tangency is called the Chapman- Jouguet point 
(p,, -,,c ) , and the explosion at this point has several 

important properties. The C-J state detonation, from some 
specified initial state (Pq/''^q)^ a minimum value for the 

propagation velocity D, relative to the unreacted material, 
and a stationary minimum value of entropy. The Hugoniot 
curve is tangent to the isentrope through the C-J point, and 
therefore the Rayleigh line is also tangent to that same 
isentrope. The flow of the reaction products relative to 
the reaction front is sonic behind a C-J detonation. The 
Chapman-Jouguet hypothesis states that the detonation wave 
for which the wave flow proceeds at a sonic velocity is the 
one that actually occurs when the reaction wave is unsupported. 
Subsequent arguments by Becker [Ref. 9] showed, through 
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entrophy considerations, that only the C-J detonation is 
stable. The C-J plane is the locus of C-J points, where the 
chemical reactions have come to thermodynamic equilibrium and 
the detonation properties predicted by the equations of state 
exist. 
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III. MODEL DEVELOPMENT 



A. GENERAL DESCRIPTION AND GEOMETRY 

This investigation is concerned with the centrally initiated 
detonation of spheres of condensed organic explosives of 
"CHNO" composition both with and without aluminvim. A standard 
mass of ten kilograms is modeled, contained in a sphere of 
about 10 to 15 centimeters radius, depending on the density 
of the explosive. The volume taken up by the initiator at the 
center is ignored as is any energy contribution from the 
initiator, since both are very small compared to those quanti- 
ties for the main charge. The explosive is detonated in a 
rigid container of irregular geometry, with an average distance 
from the charge to the wall of five meters, filled with normal 
air at standard conditions. When the detonation shock wave 
reaches the surface of the explosive, the initial shock wave 
is established in the surrounding air, there is a contraction 
of the surface of the explosive due to the reflected shock 
wave, and the reaction zone ^ the surface of the sphere. 

The shock wave displaces the air slightly away from the sur- 
face of the sphere and the compression of the surface itself 
results in a rarefaction; any chemical interaction of the 
medium with the material in the reaction zone is initially 
ignored, for reasons developed below. The reaction zone, at 
the surface, is referred to as the fireball in some contexts; 
after the hot gases produced by the detonation reaction have 
begun expansion and encounter material from the surrounding 
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medium, the surface does indeed burn in some cases. This is 
called afterburn and is most evident in oxygen-poor explosives 
or those mixed with reactive metals. For purposes of this 
model, any medixim-product reactions that occur do so after 
the beginning of the hot gas expansion and are considered to 
slow down the decay of the high temperature-high pressure 
conditions when volume expansion occurs. 

The surface of the explosive is modeled as a blackbody. 

This requires thermodynamic equilibrium, which is demonstrated 
below. Spherical geometry is selected for several reasons. 

An isotropic blackbody obeys the Stef an-Boltzmann law, and a 
spherical surface is isotropic. For most explosive configura- 
tions of interest, at some time after completion of detonation, 
the expanding mass of hot gases assiames an approximately 
spherical geometry. The true geometry of a real system of 
explosive and confining chamber can be transformed into a 
spherical coordination description. Spherical geometry 
greatly facilitates computations of radiant energies. Addi- 
tionally, the experimental results reported in the literature 
for many explosives of interest have been for spherical masses 
of explosive. 

The initiation of detonation of a spherical mass of con- 
densed high explosive from the center produces an expanding 
spherical detonation shock wave which propagates through the 
explosive. Although several sources have indicated a problem 
with maintaining a spherically diverging finite reaction zone 
of constant width [Ref. 8] , because it would of necessity 
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decrease in width as it expands, Cook [Ref. 10] employs a 
spheriform shock front model, and photographic evidence for 
such an expansion was reported by Rudlin [Ref. 11]. A point 
of consideration to resolve this quandry is that the value 
of the propagation velocity, D, builds up to the stable C-J 
value more slowly in the spherical case than in the case of 
plane wave propagation through an end-initiated cylindrical 
charge, and is a function of the radial distance out from the 
center [Ref. 8, ch. 9]. For purposes of the model used in 
this investigation, the C-J conditions are considered to exist 
at the surface of the sphere at the instant of departure of 
the initial shock wave into the medium. A minimum size exists 
below which this model does not apply, due to not reaching C-J 
conditions by the time the detonation wave has reached the 
surface, and due to an unfavorably small mass-to-surf ace area 
ratio (since heat of detonation and, therefore, the energy 
evolved, is a function of the mass of explosive and the radi- 
ant energy emitted is a function of the surface area and 
temperature) . 

B. THERM0DYNAT4IC EQUILIBRIUM 

Existence of thermodynamic equilibrium in the detonation 
reaction zone during a sensible time period is an essential 
element of the foundation of assumptions upon which internal 
blast studies have been built. Sensible time is chosen to be 
that time required for the reaction zone to move through the 
detonating explosive a distance equal to the width of the 
reaction zone . Although an exact value for the width of a 
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reaction zone is dependent on a number of variables, such as 
initial conditions of the explosive (density, packing, configu 
ration, temperature, etc.) and the chemical composition and 
experimental measurement is, at best, very difficult, values 
on the order of 0.1 millimeter have been reported, based on 
pressure distribution in the reaction zone determined by 
shockwave measurements in thin layers of inert material placed 
in contact with the end surface of the charge [Ref. 12]. In 
Ref. 10, Cook makes a correlation between the reaction zone 
width and an area of ionization which exists immediately be- 
hind the detonation front. Using this correlation, reaction 
zone widths on the order of centimeters are reported. The 
reaction zone does propagate through the explosive immediately 
behind the detonation wave, with essentially the same velocity 
D. As a typical case, TNT is considered. (Appendix A is a 
glossary of explosive acronyms; the common names of the explo- 
sives considered in this investigation are defined, and their 
compositions are given.) Reported reaction zone widths of 
2.5 cm [Ref. 10] and 0.1 mm [Ref. 12] are used to give an 

upper limit value for the sensible time of 3.7 x 10 ^s and 

— 8 

a lower limit value of 1.5 x 10 s, based on a reported 
detonation velocity value of 6830 m/s [Ref. 1] . If a large 
number of collisions occur between molecules of the product 
gases during this sensible time, it is reasonable to conclude 
that thermodynamic equilibrium exists. 

Geometric collisional diameters of product species have 
been reported in the range of one to five Angstroms (1 to 
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5 X 10 , determined by viscosity, heat conductivity and 

diffusion experiments [Ref. 13] or molecular beam apparatus 
experiments [Ref. 14] . Values of such collision diameters 
for diatomic and polyatomic molecular species are found to 
be, for example, 2.47 x lO^^^m for H 2 and 4.82 x 10”^^m for 
CO 2 / with many values reported in the range of 3.5 to 5.0 x 
10 ^^m for molecules of the species found in detonation 
product gases. Order of magnitude agreement was found through- 
out the literature examined. 

The number of particles, their collisional cross sections 
and thermal velocities can be combined to give a measure of 
the average time between collisions, through the relationship 

t =-^ 
c n a V 

where n is defined as the number density of particles in the 
gasf V as the average velocity and a as the collisional 
cross section. Reference 15 contains a form of expression 
for V (the Maxwell mean speed) ; using the variables of this 

L 

analysis, v = (SkT/m) . Reference 16 cites an equivalent 

expression, and provides a method to obtain an estimate of 

cross section for collisions between molecules of different 

2 

species if molecular radii are known: a = ir (r^ + rj^) , 

where r^ is the radius of molecule species "a” and rj^ is 
that of molecule species "b." Variations of molecular diam- 
eters, and the large number of possible varieties of colli- 
sion events that could occur between the dozen species usu- 
ally found in product gases, call for a simplifying 
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approximation. An average molecular collision diameter of 
4 X 10 ^*^m is selected, which results in a value for a of 
50 X 10 ^^m^ (5 X 10 ^^cm^) . 

Composition of product gas for a large number of organic 

chemical high explosives is reported by Mader [Ref. 1]. An 

average molecular weight for these gases was determined to 

be 30 (28.8 for TNT of density 1.60 gm/cm^ , as an example). 

The particle density, n, is determined by multiplying the 

19 3 

Loschmidt number (2.69 x 10 particles/cm ) by the detona- 
tion pressure in atmospheres and by the fraction 273/ 
(detonation temperature) . For condensed TNT of density 

3 

1.60 gm/cm , Ref. 1 reports detonation pressure of 194 kbars 

and detonation temperature of 2971 K. Using these values 

and an average molecular weight of 28.2, n was calculated to 
23 3 

be 4.73 X 10 particles/cm , v was calculated to be 1.49 

X 10^ cm/s, and a value of t was calculated to be 2.84 x 

c 

10 ^^s. Comparison of t^ to the shortest sensible time 

— 8 

value of 1.5 X 10 s demonstrates that a very large number 
of collisional events, on the order of millions, occur during 
the sensible time. Thermodynamic equilibrium thus exists in 
the reaction zone at the C-J plane. 

C . BLACKBODY APPROXIMATION 

The maximum radiant energy emitted by any isotropic body 
at thermodynamic equilibrium is, by definition, that radiated 
by a blackbody. A determination of the spectral energy given 
off by a blackbody at the detonation temperature for some 
explosive provides a measure of the maximum energy which 
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could be lost from the product gases through thermal radia- 
tion. The surface of a sphere of detonated explosive, immedi- 
ately after the departure of the shock wave into the 
surrounding medium, is modeled as a blackbody made up of hot 
gas which emits the spectrum characteristic of the tempera- 
ture of detonation (the C-J temperature) . Justification of 
this model follows order of magnitude comparisons between the 
excited state lifetimes and the time between collisions in 
the gas . 

The lifetime of an excited state is defined as the aver- 
age time an atom spends in this state before making a transi- 
tion. Radiation can be given off by the atom when transition 
to a lower state occurs. Lifetimes for many states are on 

-9 

the order of nanoseconds (10 s) [Ref. 17]. Comparison of 

-9 

typical excited state (radiator) lifetimes (10 s) to the 

average time between collisions for the gas molecules in a 
-15 

detonation (10 s) , as computed above, leads to the conclu- 
sion that the individual spectral line widths for the dozen 
or so species present in the detonation reaction zone will 
be collisionally broadened to the point of showing a continu- 
um. Consider, as an example, the case of the 10.2 micron 

13 

CO 2 line. It has a center frequency of 2.94 x 10 Hz. CO 2 

has a collision diameter of 4.84 x 10 ^^m. At a C-J tempera- 

5 

ture of 3000 K and pressure of 2 x 10 atmospheres, a time 

-15 

of collision of about 10 s is calculated. The line width 
(frequency spread) due to collisional broadening is about 
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o = 



ir T, 



[Ref. 18] . 



14 

In this example, a line width of about 3 x 10 Hz results, as 

compared to a room temperature width of about 50 MHz. This 

collisional broadening smears the 10.2 micron line from 
6 20 

about 10 Hz through 10 Hz. The line spectriam of an explosion 
is very rich [Ref. 19] ; such collisional broadening is opera- 
ting on all species line widths to some comparable degree, 
causing them to smear together into an apparent continuum. 
Within the excited state lifetime, hundreds of thousands of 
collisions occur. Thus, thermodynamic equilibrium and continu- 
um radiation are characteristic of the reaction zone, and 
application of a blackbody model is quite justified. 

The radiant exitance of a blackbody is given by the 
Stefan-Boltzmann law, 

M = a T^ where a = 1.36x10 ^^cal cm ^ K ^ s 



Using this relationship, the energy emitted by the model black- 
body can be found if the surface area, temperature and time 
of radiation are known. 

The radiant surface of the detonation is the surface of a 
sphere occupied by a ten kilogram mass of condensed explosive, 
and a function of the density of the explosive. Different 
loading techniques produce different densities for the same 
explosive. Initial density is a necessary parametric input 
to any solution of the equations of state and is reported, 
for example, by Mader [Ref. 1], The surface area for the 
reaction zone is somewhat reduced initially because of the 
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condensation which occurs in reaction to the departing shock 
wave; subsequent pressure driven expansion of the reaction zone 
will provide a compensating increase in surface area. The 
surface area of the unreacted mass is considered to be a good 
approximation. As the fireball expands, the surface area in- 
creases; the temperature, however, decreases. The rate of 
energy radiation is a linear function of the radiator area, 
but depends on the fourth power of the temperature, as seen 
above in the Stef an-Boltzmann law. Thus, as the fireball ex- 
pands, exitance decreases. Some argument could be advanced 
that the fireball would encounter additional oxygen from the 
air as it expands, thus helping to maintain the temperature 
through exothermic formation of additional oxides. This has 
been ignored for three reasons: (1) since air is only 20% 

oxygen, for every unit of mass of oxygen taken into the fire- 
ball, four units of mass of neutral (non-reacting) gas are 
also taken in and all the new material taken in must be raised 
to the equilibrium flame temperature, thereby compensating for 
any exothermicity of the oxide formation; (2) except in the 
case of very oxygen deficient or highly metallized explosives, 
formation of additional oxides may not be particularly domi- 
nant because of the decrease in pressure which would accompany 
the expansion; and, (3) the number of oxygen molecules in the 
air immediately surrounding the surface available for reaction 
is small compared to the number of product gas molecules. 
Although afterburn does indeed contribute to overall yield in 
some explosives [Ref. 20] , the consideration made is that the 
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temperature at the C-J plane (the detonation temperature) is 
the maximum, and that the exitance computed for that tempera- 
ture and for the surface area of the original explosive mass 
represents a maximum power level; total maximiam energy loss 
through thermal radiation is then established for a time 
period. The numerical value used for this time period varies 
for each explosive, dependent on available data. Actual light- 
intensity measurements are available for only some materials; 
the narrow bandwidth of the sensors used to obtain these data 
make them of marginal value. Experimental times of arrival 
of steady overpressure level or other time of significance is 
used where available. The maximum time of significance is in 
all cases determined to be the time to establish the relative 
quiescent peak overpressure which is considered the primary 
damage mechanism in the studies reported by Kinney, Sewell 
and Graham [Ref. 2] . 

D. OPACITY OF THE REACTION ZONE 

Thermal photons emitted interior to the surface do not 
contribute to the radiant flux emanating from the fireball. 

The reason for this is the opacity of the reaction zone. 

The radiation, as photons, is absorbed or scattered before it 
can cross out of the sphere through the reaction zone. Time 
between events for photon-molecule interaction is found using 
the previous relation: 

t = -i- ; 

c n a V 

V in this case is the velocity of a photon traversing the 
body of gas. As an approximation, c^, the speed of light in 
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vacuum is used. a in this case is the cross section of the 
molecules, an admittedly smaller number than the effective 
cross section seen by the radiation as an E-M wave. An order 
of magnitude comparison is again employed; a too-small a, 
which results in a too- large t^, is inconsequential to the 
outcome of the comparison since, using the average molecular 
diameter of 4 x 10 ^^m and the n from the thermodynamic equilib- 
rium calculation above, t^ is found to be on the order of 
5x10 ^^s (4.72x10 ^^s in the case of TNT of density 1.60 

3 

gm/cm , from the Mader data) . The time for the photon to cross 
the reaction zone of typical 2.0 cm width is 6.7x10 ^^s, 
minimum. Thus, it is not probable that a photon from the 
interior would not be scattered or absorbed in the reaction 
zone transit, thereby escaping to contribute to the radiant 
flux. 

E. REFLECTION AND ABSORPTION 

In this model, all thermal radiation emitted by the fire- 
ball is considered to be absorbed by container surface. 

Typical reflectivity profiles for substances used as coating 
material, such as linseed oil based paints, show a maximum 
reflectivity of about 80%, with values under 50% more common; 
reflectivity is, of course, wavelength dependent. Reference 
21 contains a number of graphs of reflectivity vs. wavelength 
for many coating materials. Radiation emitted from the fire- 
ball will encounter surfaces tens of thousands of times per 
millisecond in the container, and be partially reflected and 
partially absorbed at each encounter. After 25 reflections. 
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even wavelengths with a reflectivity as great as .8 would be 
absorbed to less than .5% of their original intensity. About 
.05% of the equatorial cross section of the container is 
composed of the fireball cross section, under the model. 
Absorption in the fireball, with subsequent re-radiation, can 
be assumed for only some small portion of the radiated energy, 
and is ignored. Atmospheric attenuation within the container 
is likewise ignored because of the small path length involved 
to virtual extinction through surface absorption. 

F. TIME OF RADIATION 

Determination of the actual length of time for which 

thermal radiation losses are a mechanism reducing peak over- 

pressvire, and at what rate, involves complexities beyond the 

scope of this investigation. As an absolute maximum, the 

amount of time, called t , that this radiant flux reduces the 

r 

internal energy available for overpressure generation is that 
from the time of arrival of the reaction zone at the surface 
of the explosive to the time of establishment of a level of 
quasi-static overpressure at the surface of the container. 
During this time, the temperature decreases the C-J tem- 

perature, although by how much is a function of the amount of 
afterburning . 

As a maximum, the total possible exitance is considered 
to be that from a blackbody, radiating at the C-J temperature, 
with surface area equal to that of the unreacted explosive 
mass, for a time, t^, as described below. Radiant power 
decreases as the fireball grows because the temperature 
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decreases with volume expansion and exitance is a function of 
the fourth power of the temperature. This power decrease due 
to temperature drop dominates the power increase attributable 
to greater radiator surface area and afterburning where 
present. 

Times of arrival of the initial shock and quasi-static 
overpressure level have been reported for TNT [Ref. 6] . These 
times are each about five milliseconds. If the shock had main- 
tained the C-J detonation velocity for the transit across the 
five meter average distance to the container wall, its time of 
arrival (a little over half a millisecond) would have been 
about one order of magnitude smaller than the measured arrival 
time. Application of the same 10:1 ratio and approximate 
equality of initial shock and quasi-static overpressure arrival 
times to other explosives is an artificiality, although not 
unreasonable for inclusion in the model since the mechanisms 
which cause velocity decay operate on both the initial shock 
and overpressure waves. When no experimental data are avail- 
able, a time of radiation is chosen to be twenty times the 
idealized initial shock transit time as computed from the 
C-J detonation velocity. Where available, actual reported 
arrival times are used. 
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IV. RESULTS 



A. DATA 

Determination of the energy of detonation contained in a 

ten kilogram mass of explosive was rather direct since values 

of q, the specific heat of detonation, are reported in the 

literature, and Q is qxm. The value of q depends, to some 

extent, on the condition of the explosive; where full data are 

available for more than one set of initial parameters for a 

given explosive, calculations were made for each set. Specific 

volume, V, the reciprocal of density, was determined, since 

the surface area is a constant times the 2/3 power of the spe- 

2/3 2 

cific volume (A = 2245 (v ) cm ) for a constant ten kilogram 

so 

explosive mass. 

As noted above, the explosive energy is partitioned, with 
some portion going into the initial shock wave. The precise 
partition is highly variable, dependent in part on the density 
of the surrounding medium. Cole [Ref. 22] reports that 53% 
of the explosive energy of TNT is transmitted to the surround- 
ing medium in an underwater explosion. Cook [Ref. 10] examines 
these variances in fluid and condensed media in some detail, 
although no empirical relations are established, per se. For 
purposes of this investigation, an examination of the ratio of 
the total radiant energy emitted to half the total explosive 
energy is deemed appropriate. Division of this (Mxt^)/Q ratio 
by the decimal fraction of Q retained in the explosive after 
departure of the initial shock wave, were it known, would 
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provide the percentage energy loss; it is not considered 
probable that this decimal fraction would be less than .50 
in air. Thus, the maximum radiant energy loss reported in 
the data represents the evaluation of (Mxt^)/.5Q. 

Data and calculation results are contained in Table I. 
Sources of the numerical values used are cited, by reference 
number in parentheses. If no reference is cited, values are 
computed from other entries in the table, using a formula or 
technique described above. The value of q for ALEX is esti- 
mated on the basis of values for similarly aluminized 
explosives . 

B. ANALYSIS 

Limited luminosity data are contained in the literature 
for comparison. Rudlin [Ref. 6, 11 and 19] reported some 
observations on radiant emissions from TNT and pentolite. 

His analysis of data from a series of multi-ton TNT explo- 
sions, and a scaled comparison to data obtained on a series 
of eight pound sphere detonations are inconclusive due to 
an absence of quantitative measurements of thermal radiative 
power for the smaller charges. In addition, the instrumenta- 
tion used on the smaller masses was band centered toward the 
shorter wavelengths in the visible and near infrared region 
of the spectrum. Rudlin observed that there was substantial 
energy in the infrared not accounted for in his investigation, 
due to the instrxamentation red cut-off, at about .7 microns. 
Some data are available for multi-ton TNT events which give 
an integrated, broad bandwidth radiant energy measurement 
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and percentage of explosive yield [Ref. 19] . Unfortunately, 

to quote Rudlin [Ref. 19], "scaling from small charges to 

multi-ton charges, does not, in general, hold." Of interest 
is that luminosity times reported are on the order of a milli- 
second for the smaller charges. Reference 17 also contains 
luminosity data. Although no units of power or energy are 
given, luminosity times in a no-afterburn environment (N 2 
filled explosion chamber) are also very short, again on the 
order of a millisecond or two for lightly aluminized or non- 
aluminized explosives. 

The times of radiation used in Table I are longer than 
those reported for the small explosive charges and may result 
in an overestimate of the total radiant energy emitted. A 
precise correlation is not possible because of the fact that 
the sensors used to determine the reported experimental times 
were not sensitive to the infrared radiation which makes up 
the majority of the radiant energy emitted. Since maximum 
attenuation effect determination is the object of this study, 
the t^ is, within the other model criteria, considered to be 
of an appropriate magnitude. 



33 



Table i 



Explosive 


( cm^ ) 


<1 

(cal/gm) 


D 

(m/s) 


rK ) 


Loss 


(cm3/gm) 


Q 

(Meal) 


(ms) 


k 

Itecal/s ) 


(/O 


RBX 


1.78 (7) 


1480 (7) 


8000 (1) 


2970 (1) 


.03 


• 562 


14.8 


12.2 (23) 


0.162 


PSTN 


1.67 (1) 


1510 (7) 


8180 (10) 


3400 (1) 


.05 


.578 


15.1 


12.2 (10) 


0.283 


TNT 


1.54 (7) 


1290 (7) 


6950 (1) 


2940 (1) 


.03 


.649 


12i9 


11.3 (6) 


0.171 


HM2 


1.89 (7) 


1480 (7) 


9124 (23) 


2364 (1) 


.01 


.529 


14.8 


11.0 (23) 


0.0624 


COMP B 


1.69 (7) 


1400 (7) 


7840 (23) 


2770 (1) 


.02 


.592 


14.0 


12.8 (23) 


0.127 


Tetryl 


1.71 (7) 


1450 (7) 


7850 (23) 


2917 (1) 


.03 

1 


.585 


14.5 


12.7 (23) 


0.156 


TATNB 


1.89 (1) 


1080 (7) 


7500 (23) 


2128 (1) 


.01 


.529 


10.8 


13.3 (23) 


0.041 
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T/uBLS I 
( cont . ) 



i..annitol 


1.73 (23) 


1390 (10) 


1 

8260 


(23) 


6040 (10) 


.49 


Hexanitrate 


.578 


13.9 


12.1 


(23) 


2.82 


DATNB 


1.80 (7) 


1150 (7) 


7559 


(1) 


2667 (1) 


.02 


.555 


11.5 


13.2 


(1) 


0.104 


AI^ 20 


1.80 (1) 


1800 


7451 


(1) 


5142 (1) 


.21 


.556 


18.0 


13.4 


(1) 


1.44 


Tritonal 


1.42 (10) 


1100 (10) 


4930 


(10) 


4100 (10) 


.25 


80-20 


.704 


11.0 


20.3 


(10) 


0.683 


Tritonal 


1.79 (10) 


1750 (10) 


7020 


(10) 


5410 (10) 


.29 


80-20 


.559 


17.5 


14.2 


(10) 


1.77 


45/30/25 

RBX-TUT- 

Aluminimi 


1.51 (10) 


1060 (10) 


5470 


(10) 


4400 (10) 


.29 


.662 


10.6 


17.4 


(10) 


0.86^ 


45/30/25 

* RBJC-TNT- 
1 

j Aliominm 


1.82 (10) 


1630 (10) 


6760 


(10) 


5700 (lO) 


.3^ 


.549 


16.3 


14.8 


(10) 


2.16 
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V. CONCLUSIONS 



The basic assumption made in previous investigations of 
internal blast [Ref. 2 , 3, 4 and 5] that explosion products 
are in thermal equilibrium is shown to be valid. By utiliza- 
tion of a blackbody radiator model/ it is shown that explosive 
energy losses through thermal radiation from the surface of 
the explosion fireball are small, ranging from about .5% to 
under .01%. Demonstration of the validity of a second major 
assumption made in these prior investigations, that thermal 
radiative losses could be ignored, is thereby achieved. 

With the recent development of detector systems which 
operate in the infrared, better experimental data on the 
spectral intensity profile of detonating explosives may be 
obtainable to the end of experimentally verifying detonation 
temperatures predicted by thermo-hydrodynamic codes. It is 
probable, however, that the experimental error induced by the 
apparent irreproducibility of explosive results to within a 
tolerance finer than five percent shall remain a major impedi- 
ment to consistently accurate prediction of explosive 
parameters . 
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APPENDIX A 





GLOSSARY OF EXPLOSIVE ACRONYMS 


ALEX 20 


Aliiminized Explosive; 20% alxominum in 80% 
Composition B; empirical formula, based on 
100 grams per mole, is 

^1. 873^2. 469^1. 613°2. 039^^7. 338* 


Comp B 


Composition B; mixture of 64% RDX and 36% 
TNT; empirical formula, based on 347.07 
grams per mole, is 

*^6.85^8, 75^7. 65°9. 30* 


DATNB 


Diamino trinitrobenzene; chemical formula is 

C.H(.NcO,.. 

6 5 5 6 


HMX 


Cyclotetramethylene tetranitramine; chemical 
formula is C^HgNgOg. 


PETN 


Penta erythritol tetranitrate; chemical 

formula ia CcH^N.O, 

5 6 4 12 


RDX 


Cyclotrimethylene trinitramine; chemical 
formula is C,H,N-0^. 


TATNB 


Triamino trinitrobenzene; chemical formula 

is C^H-N^O^. 

6 6 6 6 


Tetryl 


Trinitrophenyl methyl nitr amine; chemical 
formula is C-H-NpOo. 

/ D D O 
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TNT 



Tritonal 



Trinitrotoluene; chemical formula is 
C 7 H 5 N 30 g. 

Aluminized TNT; numbers index percentages 
of TNT and aluminum (e.g., 80-20 is 80% 
TNT with 20% alviminum) . 
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